Noise, coherent fluctuations, and the onset of order in an oscillated granular fluid 
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We study fluctuations in a vertically oscillated layer of grains below the critical acceleration for 
the onset of ordered standing waves. As onset is approached, transient disordered waves with a 
characteristic length scale emerge and increase in power and coherence. The scaling behavior and 
the shift in the onset of order agrees with the Swift-Hohenberg theory for convection in fluids. 
However, the noise in the granular system is four orders of magnitude larger than the thermal noise 
in a convecting fluid, making the effect of granular noise observable even 20% below the onset of 
long range order. 



Recent experiments Q and simulations Q on flows of 
dissipative granular particles have been found to be de- 
scribed by hydrodynamic theory, although the granular 
systems exhibit much larger fluctuations than fluids - a 
single wavelength in a pattern in a vibrated granular layer 
might contain only 10 2 particles instead of the 10 22 par- 
ticles in one wavelength of a pattern in a vibrated liquid 
layer. Thus fluctuations play a more significant role in 
granular systems than in fluid flows. In fluids the fluctua- 
tions are driven by thermal noise and are described by the 
addition of terms to the Navier-Stokes equations 0, • 
This fluctuating hydrodynamic theory has been found 
to describe accurately the dynamics near the onset of a 
convection pattern in a fluid 0, 0, • Here we describe 
an experimental study of fluctuations near the onset of 
square patterns in granular layers, and we demonstrate 
the applicability of fluctuating hydrodynamics to this in- 
herently noisy system. 

Experiment — We study a vertically oscillating layer 
of 170 /im stainless steel particles as a function of T, the 
peak plate acceleration relative to gravity; the container 
oscillation frequency fd is fixed at 30 Hz ||. The layer 
(depth 5 particle diameters) fluidizes at T w 2 0, and 
a square wave pattern with long range order emerges at 
r w 2.77. The layer is illuminated at a low angle and 
fluctuations in the surface density are observed as fluctu- 
ations in the light intensity recorded on a 256x256 
pixel CCD camera. 

Coherent fluctuations — Fluctuations are evident in 
the snapshot of the layer shown in Fig. ^a). Spatial 
coherence with no orientational order is revealed by the 
ring in the spatial power spectrum S(k x ,k y ) shown in 
Fig. ^b). The increase in the spatial and temporal co- 
herence and the power of the fluctuations with increasing 
r is illustrated in Fig. [21 where insets with each spectrum 
S(k x ,k v ) show the corresponding azimuthally averaged 

structure factor, S(k) — (S(k Xl k y )) g with k — ^jk\ + ky. 

(There is slight (< 1%) hysteresis in all measured quan- 
tities between increasing and decreasing T, but we will 
only discuss increasing T.) The power of the dominant 
mode increases while its width decreases with increasing 
r [llj. The noise is readily visible at T = 2.2, which is 
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FIG. 1: (a) Snapshot of an area 6.25 x 6.25 cm 2 in a con- 
tainer oscillating with F = 2.6. (b) The spatial power spec- 
trum of (a) has an intense ring corresponding to randomly 
oriented spatial structures with a length scale of 0.52 cm (for 
improved signal-to-noise, 100 spectra were averaged to obtain 
the spectrum shown), (c) Space-time diagram for the row of 
pixels in the box in (a); the period of the localized transient 
oscillations is 2T = 2/ fd- 



20% below the onset of long range order, while in high 
sensitivity experiments on Rayleigh-Benard convection, 
the noise became measurable only at 1% below the onset 
of convection |]| ■ 

Local transient waves — Transient localized structures 
oscillating at /d/2 are visible in Fig. dc). Temporal 
power spectra of the intensity time series for each pixel 
in the images are averaged to obtain the power spectra 
S(f) shown in Fig. El The peak at /d/2 emerges for T > 2 
and increases in power with increasing T. The half-width 
at half-maximum of the peak, denoted 1/r, decreases as 
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FIG. 2: Spatial power spectra S(k x , k y ) and temporal spectra 
S(f) and structure factors S(k) (insets) show that the fluctu- 
ations increase in power and coherence as F increases from 2.2 
to 2.6, and at T — 2.8, above onset, a coherent square spatial 
pattern has emerged. S(k x ,k y ), S(k), and S(f) have been 
divided by 10 at V — 2.8, where there is a strong peak at 15 
Hz corresponding to the spatial pattern (this peak mixes with 
the 30 Hz drive frequency to give a peak at 45 Hz). The value 
of k corresponding to a maximum in the structure function 
for r = 2.8 corresponds to a length of 0.52 cm. 



r increases, indicating increasing temporal coherence of 
the noisy structures. Except for the peak at fd/2, S(f) 
does not depend strongly on T for T > 2. The shape of 
S(f) is similar to that for surface excitations (ripplons) 
on a fluid driven by thermal noise [l2|. 

Fluctuating hydrodynamics — The phenomena we have 



described have the features of noise-driven damped hy- 
drodynamic modes close to a bifurcation. We now in- 
terpret the observations using the Swift-Hohcnbcrg (SH) 
model, which is based on the Navier-Stokes equation and 
was developed to describe noise near the onset of long 
range order in Rayleigh-Benard convection 0. Swift- 
Hohenberg theory predicts that below the onset of or- 
dered patterns, noise drives a ring of modes that in- 
creases in power as onset is approached; our observa- 
tions in Fig. [21 are in qualitative accord with this pre- 
diction. The theory also predicts that the nonlinearity 
of the fluid acting on the noise will lead to an increase 
of the critical value of the bifurcation parameter for the 
onset of long range order. The observations also agree 
qualitatively with this prediction: for the stainless steel 
particles the patterns are noisier than those obtained in 
previous experiments on lead particles, which are more 
dissipative; further, the onset of long range order for the 
stainless steel particles occurs for T « 2.77, while for 
lead, r«2.5@. 

The SH model describes the evolution of a spatial 
scalar field ip(x,t), 

^ = (e-(^ + ktr)^-^ + n(x,t), (1) 

where e is the bifurcation parameter and rj is a stochastic 
noise term such that (?y(x, t)r](x',t')) = 2F5(x — x')S(t — 
t') 1 where F denotes the strength of the noise. In the 
absence of noise [F = 0), called the mean field (MF) ap- 
proximation 0, , there is an onset of stripe patterns 
with long-range order at e = e" F = 0. (Our experi- 
ments yield squares at pattern onset with slight hystere- 
sis, but we compare our observations to the the simplest 
model for noise interacting with a bifurcation, Eq. ^ 
which yields stripes at onset via a forward bifurcation; 
a more complicated model yielding square patterns and 
hysteresis is described in |l5|.) For F ^ 0, the onset 
of long-range (LR) order is delayed until e = eJ; R > 0, 
where eJ: R oc F 2 / 3 . For < e < the pattern is disor- 
dered, and appears cellular-like |l6(. We define the delay 
in onset as Ae c = e^ R - e" F . 

To compare results from the experiments with re- 
sults obtained by integration of the SH model, we must 
first compute the reduced control parameter e = (r — 
rjf F )/r" F . However, we have no a priori way to deter- 
mine r" F since the theory predicts a smooth change in all 
quantities as the mean field bifurcation is crossed. We de- 
termine T" F by performing a two-parameter fit to obtain 
agreement between experiment and the SH model for the 
maximum value of S(k) for all e below onset and for the 
critical value for the emergence of long range order e^f as 
determined from the onset of angular correlation in the 
radially averaged structure factor (C(9 = ir/2) for the 
emergence of squares in the experiments and C{9 = tt) 
for stripes in the SH model) . In the fitting procedure we 
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FIG. 3: The growth of the noise power and the onset of 
long-range order in (a) experiment and (b) Swift-Hohenberg 
model. The log of the maximum of S(k) (o) increases through 
the mean field onset (e" F = 0), while the onset of long range 
order, indicated by appearance of angular correlations of the 
radially averaged structure factor [C(9 = it/2) for the ex- 
periment and C(9 = n) for SH equation (•)], is delayed to 
£c R ~ 0.04. The integration of Eq.^uses a scheme described 
in p^| : the solution is obtained on a 128 x 128 grid with fco = 1 
and integration time-step 0.5. 



vary Tf F for the experimental data and F in Eq. ^ .18]. 
Thus, in addition to a prediction of T" F , the procedure 
allows us to estimate the size of the noise required to shift 
e™ F = up to the observed eJ; R , i.e., to shift r™ F up to 

r^ R = 2.77. 

The two-parameter fit yields T" F = 2.64 and F = 
3.5 x 10 -3 , giving a delay in onset of Ae c = 0.04 (Fig. 01. 
The noise strength is a factor of 10 4 larger than the noise 
in the convecting fluid experiments of Wu et al. [f| and is 
a factor of 10 larger than in recent experiments on con- 
vection near a critical point [l9| . The theory also predicts 
that at eJ; R there should be a jump in (ip 2 ) proportional 
to F 2 / 3 ; however, the predicted jump, only 10~ 3 in (ip 2 ), 
is too small to detect with the precision of our exper- 
imental measurements. (Also, such an effect would be 
dominated by the hysteresis in the transition to square 
patterns.) Now that we know the appropriate reduced 
control parameter e for the experiment, we can test pre- 
dictions of scaling given by Eq. ^ 

Scaling — Linear theory predicts the following scal- 
ings for thermal convection with small noise levels near 
the onset of long range order: both the noise peak inten- 
sity (at the wavenumber selected by the system) and the 



correlation time should scale as |e| 1 and the power in 
the fluctuations should scale as |e| - ^ 0j@- The scaling 
behavior found in both the experiment and SH model for 
large noise levels is shown in Fig. 01 where the noise peak 
intensity was determined from the maximum of S(k) (see 
insets of Fig. EJ) , and the power in the fluctuations was 
determined from the area under the peak in S(k), af- 
ter subtracting the approximately constant background 
due presumably to incoherent grain noise. The Swift- 
Hohenberg theory agrees remarkably well with the ob- 
servations. As expected, both experiment and theory de- 
viate considerably from the linear theory prediction for 
small |e|, where nonlinear effects are large, but for large 
|e| the experiment and simulation approach the scaling 
predicted by the linear theory. Finally, we have deter- 
mined from the experimental data the correlation time 
for the patterns by fitting the fd/2 peak of S(f) to a 
Lorentzian and computing the half-width at half maxi- 
mum, 1/r. Again, for large e the observations approach 
the scaling predicted by linear theory, but deviate from 
linear theory for small e (Fig.0J. 
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FIG. 4: The noise peak intensity, total noise power, and the 
noise coherence time as a function of the distance e below 
the onset of patterns. The measured quantities (•) compare 
well to the Swift-Hohenberg equation (solid curve) with F = 
3.5 x 10 -3 . Results from both experiment and SH theory 
approach the predictions from the linear theory (shown by 
the dashed lines with the values for their slopes) when e is 
large, far below the onset of long range order. 

Conclusions — We have shown that a vertically oscil- 
lated layer of grains exhibits behavior consistent with the 
theory of fluctuating hydrodynamics for Navier-Stokes 
fluids. This indicates that fluctuations must be included 
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in the hydrodynamic equations for granular media [20j- 
In fact, the large fluctuations present in granular fluids 
suggests that such systems could be useful to study the 
effects of noise in nonequilibrium fluids far below a bifur- 
cation |2llE|. 

Like numerical studies of elastic gases j2^| , our exper- 
iments and simulations show that fluctuating hydrody- 
namics can apply down to length scales of only a few 
mean free paths. Fluctuations are important in fluids at 
the nanoscale, which are of current interest [241 15^ . The 
fluctuations are difficult to study in gases and liquids 
but can be studied easily in granular materials, which 
may demonstrate some essential features of the nanoscale 
flows. 
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